Thermal expansion in small metal clusters and its impact on the electric polarizability 
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The thermal expansion coefficients of NaN clusters with 
8 < N < 40 and AI7, Al]~ 3 and Al]~ 4 are obtained from ab ini- 
tio Born-Oppenheimer LDA molecular dynamics. Thermal 
expansion of small metal clusters is considerably larger than 
that in the bulk and size-dependent. We demonstrate that 
the average static electric dipole polarizability of Na clusters 
depends linearly on the mean interatomic distance and only 
to a minor extent on the detailed ionic configuration when the 
overall shape of the electron density is enforced by electronic 
shell effects. The polarizability is thus a sensitive indicator 
for thermal expansion. We show that taking this effect into 
account brings theoretical and experimental polarizabilities 
into quantitative agreement. 

PACS: 36.40.Cg, 65.70.+y, 33.15.Kr 



Since electronic shell effects were put into evidence in 
small metallic systems metal clusters have conti- 

nously attracted great interest both experimentally and 
theoretically [p]-p~0[| . Besides technological prospects, one 
of the driving forces for this research has been the funda- 
mental question of how matter develops from the atom to 
systems of increasing size, and how properties change in 
the course of this growing process. In some cases it has 
been possible to extract detailed information from ex- 
periments done at low temperatures JllJ and the related 
theories [ p"2| . In many cases, however, a deeper under- 
standing is complicated by the finite temperature which 
is present in most experiments due to the cluster pro- 
duction process, see e.g. the discussion in |f]. Whereas 
a lot of theoretical information about finite temperature 
effects in nonmetallic systems has been gained in the last 
years [|l4|], only little is known about it in metallic clus- 
ters. Here, sodium is a particularly interesting reference 
system because of its textbook metallic properties and 
the fact that it has been extensively studied within the 
jellium model, see e.g. fl5|| for an overview. Aluminum, 
on the other hand, is of considerable technological inter- 
est. Some advances to study temperature effects in metal 
clusters including the ionic degrees of freedom were made 
using phenomenological molecular dynamics (ilf | , a tight - 
binding hamiltonian | p7[ , the Thomas-Fermi approxima- 
tion [Q or the Car-Parrinello method |l^] . Recently, it 
has also become possible to study sodium clusters of con- 
siderable size p0| using ab initio Born-Oppenheimer, lo- 
cal spin density molecular dynamics (BO-LSD-MD) pi| . 

In this work we report on the size dependence of a 
thermal property which is well known for bulk systems, 



namely the linear thermal expansion coefficient 

■ 1 dT 



(i) 



For crystalline sodium at room temperature, it takes |2^] 
the value 71 x lCT 6 /-^ 1 , for Al 23.6 x lCT 6 /^ 1 . To 
the present date, however, it has not been known how 
small systems are affected by thermal expansion. At first 
sight, it is not even obvious how thermal expansion can 
be defined in small clusters. Whereas in the bulk it is no 
problem to define the length I appearing in Eq. (Jlj), e.g. 
the lattice constant, it is less straightforward to choose a 
meaningful I in the case where many different ionic ge- 
ometries must be compared to one another. For small 
metal clusters, the latter situation arises because of the 
many different isomers which appear at elevated temper- 
atures. 

We have calculated the thermal expansion coefficients 
for Nas, Naio, Naj.2, Nai4, Na2o an d Na4o in BO-LSD- 
MD simulations. Results concerning isomerization pro- 
cesses in these simulations have been presented in [E3|, 
and the BO-LSD-MD method is described in detail in 
Ref. [pi]. A meaningful length to be used in Eq. ([!]) if it 
is applied to finite systems with similar overall deforma- 
tion is the mean interatomic distance 
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where R^ are the positions of the N atoms in the clus- 
ter. Obviously, Z m i a d measures the average "extension" 
of a clusters ionic structure, and we calculated it for all 
configurations obtained in a BO-LSD-MD run. Two dif- 
ferent methods were used to calculate (3. First, we discuss 
the heating runs, in which the clusters were thermalized 
to a starting temperature and then heated linearly with 
a heating rate of 5K/ps and a time step of 5.2 fms. l m i a d 
was recorded after each time step. In this way, for Nas 
the temperature range from about 50 K to 670 K was 
covered, corresponding to 24140 configurations, for Naio 
from ca. 150 K to 390 K (9260 configurations), for Nai4 
from ca. 50 K to 490 K (17020 configurations), for Na 20 
from ca. 170 K to 380 K (8000 configurations), and for 
Na 40 from ca. 200 K to 400 K (7770 configurations). 

Fig. [l] shows how Z m i a d changes with temperature for 
Nas and Naio- Both curves show large fluctuations, as 
is to be expected for such small systems. However, one 
clearly sees a linear rise as the general trend. We there- 
fore made linear fits to the data for each cluster in two 
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FIG. 1. Mean interatomic distance in ag versus tempera- 
ture in K for Nag and Naio. The dashed lines indicate linear 
fits to the complete set of data, see text for discussion. Note 
the different slopes for the two clusters. 



ways. The first column in the left half of table | gives the 
linear thermal expansion coefficients which we obtained 
from fitting the data in the temperature interval between 
200 K and 350 K, i.e. around room temperature, where 
bulk sodium is usually studied. In order to allow for an 
estimate of the statistical quality of the fits in view of 
the fluctuations, the second and third column in the left 
half of Tabic | list the ratio of the fit parameters, i.e. 
the axis interception a and the slope b, to their stan- 
dard deviations. It becomes clear from these results that 
thermal expansion in the small clusters is considerably 
larger than that in the bulk. This can be understood as 
an effect of the increased surface to volume ratio in the 
finite systems. However, the expansion coefficient also 
strongly depends on the cluster size. This can even be 
seen directly from the different slopes in Fig. |l|. As we 
will show below, this size dependence has far reaching 
consequences for the interpretation of experimental data 
which is usually measured on hot clusters, as e.g. the 
static electric polarizability. 

In addition to the values given in Table |], we calcu- 
lated the expansion coefficient of Nai2 with a different 
method. In two separate runs, the cluster was thermal- 
ized to temperatures of about 200 K and 350 K, and 
then BO-LSD-MD was performed for 5 ps at each tem- 
perature, i.e. without heating. From the average i m i a d 
found in the two simulations, /3Nai 2 = 2.5/?buik was cal- 
culated. Thus, also the second method leads to a p that 
is larger than that of the bulk, i.e. it confirms the results 
of the heating runs. 

The average thermal expansion coefficient for the full 



TABLE I. Left half, first column: Linear thermal expan- 
sion coefficient of small Na clusters in the temperature inter- 
val between 200 and 350 K, given in terms of the bulk value 
71 x lO -6 ^ -1 . Columns two and three give the ratio of the 
axis interception a and the slope b to their standard deviations 
as obtained from the fits. Right half: Expansion coefficient 
averaged over 50-670 K for Na 8 , 150-390 K for Naio, 50-490 
K for Nai 4 , 150-460 K for Na 20 , and 200-300 K for Na 40 . See 
text for discussion. 



temperature range covered in each simulation is obtained 
from a fit to the complete set of data, shown as a dashed 
line in Fig. |l] for Nas and Naio ■ This average is of inter- 
est because it covers several hundred K for each cluster 
in the range of temperatures which are to be expected 
for clusters coming from the usual supersonic expansion 
sources |2j]. The right half of table | lists these average 
expansion coefficients and their statistical deviations in 
the same way as before. As is to be expected, the val- 
ues differ from the previous ones for the small clusters, 
because the expansion coefficient is influenced by which 
isomers are or become accessible at a particular tempera- 
ture, i.e. especially at low temperatures it is temperature 
dependent. In Fig. I] one e.g. sees from comparison with 
the average dashed line that for temperatures between 50 
K and 100 K, the thermal expansion is smaller than that 
seen for higher temperatures. However, once the cluster 
has reached a temperature where it easily changes from 
one isomer to another, the thermal expansion coefficient 
becomes nearly independent of the temperature. In the 
case of Nas, e.g., P changes only by about 5 % in the 
interval between 300 K and 670 K. 

Detailed previous investigations pp| , p3| have shown 
that small clusters do not show a distinct melting tran- 
sition. However, the largest cluster studied here, Na4o, 
shows a phase transition above 300 K [M. At the melt- 
ing point, the octupole and hexadecupole deformation of 
the electronic density sharply increase. If Z m i a d is a rele- 
vant indicator for structural changes, then melting should 
also be detectable from it. Indeed we find a noticeable 
increase in / m i ac i at 300 K, and similar fluctuation pat- 
terns as in the multipole moments. In our simulation, 
we could only determine the expansion coefficient for the 
solid phase, and it is given in the right half of table [jj 

As seen in Fig. [I], Nag shows thermal expansion already 
at 50 K. This raises the question at which temperature 
the expansion actually starts, i.e. where anharmonic ef- 
fects in the ionic oscillations will start to become impor- 
tant. In this context we note that one can compare the 
^miad at T=0 K found by extrapolation from the heating 
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data to the Z m i a d which is actually found for the ground 
state structure at T=0 K. We have done this for Nas, 
Naio an d Nai4, where the ground state structures are 
well established. In all cases, the differences between the 
two values were less than 1%. This indicates that the an- 
harmonic effects for Na clusters are important down to 
very low temperatures. Furthermore, the anharmonici- 
ties should also be observable in the heat capacities [^o) , 
where they will lead to deviations from Dulong-Pctit's 
law. We have checked this and indeed found deviations 
between 8 % (Na 20 ) and 19 % (Na 8 ) from the Dulong- 
Pctit value. 

As an example for the considerable influence of ther- 
mal expansion on measurable physical properties we dis- 
cuss the average static electric dipole polarizability a, 
which is defined as the trace of the polarizability ten- 
sor. It was one of the first observables from which the 
existence of electronic shell effects in metal clusters was 
deduced |IJ , and it has been measured for clusters of var- 
ious sizes and materials jl(J. For Na clusters with up 
to eight atoms, the polarizability was also calculated in 
different approaches [pH^ppCjl ■ These calculations qual- 
itatively reproduce the experimentaly observed trends, 
but they all underestimate the measured value. We show 
that this discrepancy is to a large part due to the fact 
that the calculations were done for T=0, whereas the 
measurement is done on clusters having temperatures of 
about 400 to 600 K ||. 

For various, different isomers obtained in our heating 
runs for Nas and Naio, we have calculated the polariz- 
ability from the derivative of the induced dipole moment 
with respect to the electric field (finite field method). 
Since highly unsymmetric isomers from the high temper- 
ature part of the simulations were taken into account, 
the full tensor was computed by numerically applying 
the dipole field in the different directions in seperate cal- 
culations. We have checked that the used field strength 
of 5 x 10~ 5 e/aQ is large enough to give a numerically 
stable signal and small enough to be in the regime of 
linear response. In Fig. || we have plotted the thus ob- 
tained polarizabilities versus Z m iad) and show three in- 
stances of ionic geometries for each cluster that demon- 
strate how different the structures actually are. Never- 
theless, within a few percent the polarizabilities are on a 
straight line. This shows that the average polarizability 
depends mainly and strongly on the mean interatomic 
distance, and only to a minor extent on details in the 
ionic configurations. Of course, the situation might be 
more complicated for clusters where the overal shape, 
i.e. the lowest terms in the multipole expansion of the 
valence electron density, is not stabilized by electronic 
shell effects. For the present clusters, however, the de- 
formation induced by the electronic shell effects persists 
even at elevated temperatures. That a is less sensitive 
to the detailed ionic configuration than, e.g., the pho- 
toabsorption spectrum, is understandable because it is 
an average quantity. 

The dependence of the polarizability on the mean 
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FIG. 2. Static electric dipole polarizability versus mean in- 
teratomic distance for different isomers of Nag (upper) and 
Naio (lower). Three examples of different geometries are 
shown as insets for both sizes. 



interatomic distance has the consequence that a also 
strongly depends on the temperature. From Fig. |^ one 
deduces that an average bondlength increase of 1 ag in 
Nas and Naio leads to an increase in the polarizability 
of about 25 A 3 . Thus, neglection of the thermal expan- 
sion in T=0 calculations leads to polarizabilities which 
are smaller than the ones measured on clusters coming 
from supersonic expansion sources . Of course, also 
underestimations of the cluster bond lengths that are due 
to other reasons will directly appear in the polarizabil- 
ity. With the Troullier-Martins pseudopotential, e.g. the 
BO-LSD-MD underestimates the dimer bond length by 
4.5%, and it is to be expected that the situation is similar 
for the bond lengths of larger clusters. Taking this into 
account, one can proceed to calculate the polarizability 
for clusters with a temperature corresponding to the ex- 
perimental one of about 500 K . In the experiments 
the clusters are spending about 10 _4 s in the deflecting 
field from which the polarizability is deduced, i.e. the ex- 
perimental timescale is orders of magnitude larger than 
the timescale of the fluctuations in the mean interatomic 



3 



distance (see Fig. [|). Thus, the fluctuations will be av- 
eraged over and can be neglected. From the average ex- 
pansion coefficients we obtain a bond length increase of 
0.48 a for Na 8 and 0.87 a for Nai at 500 K, which in 
turn leads to an increase in the polarizability of 12 A 3 
and 23 A 3 , respectively. The resulting polarizabilities of 
130 A 3 for Na 8 and 172 A 3 for Naio compare favourably 
with the experimental values 134±16A 3 and 190±20A 3 
10 



For all other cluster sizes, the two experiments 
[1,10 give different values for the polarizability. From the 
present work it becomes clear that differences in the ex- 
perimental temperatures might be the reason for the dis- 
crepancies. Therefore, an accurate measurement of the 
clusters' temperatures is necessary before further quan- 
titative comparisons can be made. However, a detailed 
comparison to both experiments showed that the theo- 
retical T=0 polarizability of all isomers underestimates 
both experimental results [ p5[ . Thus, the increase in a 
that is brought about by thermal expansion will lead to 
better agreement between theory and experiment for all 
cluster sizes. 



Thermal expansion is also observed in aluminum clus- 
ters. For AI7 we performed 5 ps of BO-LSD-MD at each 
of the fixed temperatures 100 K, 300 K, 500 K and 600 
K, for A17/ 3 at 260 K, 570 K and 930 K, and for Al7/ 4 at 
200 K, 570 K and 900 K, in analogy to the procedure 
for Nai2- From the average l m iad at each temperature, 
we calculated the expansion coefficients /3ai 7 = l-3/3t>uik, 
/3 41 - = 1.4/3buik, /?ai- = l-4/3buik- It should be noted 

Al 13 rtl 14 

that with AV[ 3 we have chosen an electronically as well as 
geometrically magic cluster pq , i.e. a particularly rigid 
one, and the fact that it also shows a larger expansion 
coefficient than the bulk is further evidence for the con- 
clusion that the increased expansion coefficient is indeed 
a finite size effects. A noteworthy difference between Al 
and Na is seen in the temperatures where the expansion 
sets in. Whereas for Na this temperature is below 50 K, 
we observe that AlJ~ 3 and AlJ~ 4 show no expansion below 
300 K. 

In summary, we have calculated thermal expansion co- 
efficients for small metal cluster and demonstrated that 
thermal expansion in these systems is larger than that 
in the bulk. For the case of sodium, the dependence of 
the expansion coefficient is not monotonous according to 
the cluster size. We showed that the average static elec- 
tric dipole polarizability of clusters whose overall shape 
is fixed by electronic shell effects depends linearly on the 
mean interatomic distance. Thus, thermal expansion in- 
creases the static electric polarizability, and we demon- 
strated that taking this effect into account brings the 
theoretical values in a close agreement with the experi- 
mental ones. 
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